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Abstract 

The quantum cosmological version of a nonsingular Universe presented by Mukhanov and Brandenberger 
in the early nineties has been developed and the Hamilton Jacobi equation has been found under semiclassical 
, (WKB) approximation. It has been pointed out that, parameterization of classical trajectories with semiclas- 

sical time parameter, for such a classically constrained system, is a nontrivial task and requires Lagrangian 
formulation rather than the Hamiltonian formalism. 
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^ • 1 Introduction 

'— ' i It transpires from Hawking-Penrose energy condition that the singularities of the general theory of relativity are 
^ unavoidable classically. In the absence of a complete and satisfactory theory of quantum gravity, it is not clear what 
^ would be the nature or the fate of singularities in the quantum domain. However, it might be possible in principle 
fS| to get rid of the singularities classically by imposing some quantum mechanical bounds on certain quantities. At 
, the singularity, some of the curvatures, eg., energy momentum and even the Riemann tensors diverge. Further, 
' quantum field theory predicts that, the scattering cross sections become infinity when all radiative corrections 
I are taken into account. Renormalization makes individual terms manageable, but the entire series diverges. So 
in order to get rid of the singularities, curvature invariant terms must be constrained to take some upper limits. 
If Planck's length Ipi is assumed to be the fundamental length, below which no length is measurable, then from 
dimensional argument, curvature invariants must be bounded as, \R\ < l^j^ , \RpiuR^'^\ < , etc. However, 
curvature invariants are infinite in number and there is no guarantee that imposing constraint on some lower order 
curvature invariant terms, all the higher order terms will be bounded. Hence, 'Limiting Curvature Hypothesis' 
(LCH) [1] had been invoked. LCH states that a finite number of curvature invariants should be bounded and one 
^ ' of those, say I2 , should have the property that /2 = singles out a particular nonsingular solution as the only 
one. In the process all the infinite set of curvature invariants are automatically bounded. If the limiting space is 
de-Sitter, then for a isotropic Universe, LCH implies that the initial and the final stages of cosmological evolution 
for a closed Universe will give de-Sitter phase, while baby Universes will arise at the centre of the Black-Hole and 
eventually gravitational collapse is avoided. 

Since Weyl tensor vanishes for an isotropic Universe, so in its simplest form, LCH can be realized by constraining 
only a couple of lower order curvature invariants, viz., and h = y/ARf.^Rt^" - R"^ . Thus in its simplest 

form, ie., for homogeneous and isotropic space-time, Mukhanov and Brandenberger [2] followed by Brandenberger, 
Mukhanov and Sornborger [3] had constructed the following action (1), which preserves LCH. 
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S = + <Pi)R - {<P2 + V34>l)^Jm^..R^•'-R^) + VMi) + V2{cf>2)]V^ d^x. 

Here, the scalar field 02 has been introduced to obtain nonsingular solutions, while (pi bounds the curvature 
invariants for these nonsingular solutions. The potentials Vi{(f)i) and V2{4>2) should be chosen in a manner such 
that (i) for 0^ <C 1 , the leading terms in the action gives back Einstein's theory and (m) their asymptotic be- 
haviour should be such that LCH may be realized, which requires — > constant as (/>i — > 00 , and V2{4>2) 
as (t>2 ^ 00. For a flat Robertson- Walker metric, they ([2] & [3]) had obtained a first order differential equation 
containing fields and the potentials and had drawn a phase diagram in (pi — 02 plane, for some particular choice 
of the potentials. The phase diagram indicates nonsingular solutions in all regions and de-Sitter phase in the 
asymptotic region. In a nut-shell, the results obtained [3] are the following. Firstly, all the homogeneous and 



isotropic cosmological solutions are nonsingular. Next, Two dimensional cosmological and Black-Hole solutions 
are nonsingular. Finally, there are evidence that nonsingular solutions may exist for four-dimensional Black-Holes 
and homogeneous-anisotropic cosmological models. Further, LCH has also been applied to dilaton-cosmology [5] 
and as a result, a class of spacially flat bouncing cosmological solutions have emerged. 

This construction is definitely appealing in itself, however, it also incorporates some additional features in quan- 
tum domain. In the path integral formulation of quantum gravity, one has to analytically continue the Lorentzian 
space of indefinite metric to the Euclidean space of positive definite metric in order to ensure convergence of the 
path integral. Nevertheless, since Euclidean action of the gravitational field is not bounded from below, the path 
integral never really converges. The action presented by Mukhanov and Brandenberger [2] and Brandenberger, 
Mukhanov and Sornborger [3J is bounded from below, and therefore, the corresponding path integral converges. 
In addition, since the initial and the final stages of cosmological evolution are de-Sitter, so one gets Lorentzian 
wormholes instead of Euclidean. Finally, as the authors had claimed that the theory can incorporate matter field 
and also can be extended to include anisotropic models, so it appears that the model is of greater interest with a 
wide spectrum of its applicability. 

In the theory under consideration, two scalar fields, one of which is non-dynamical, have been nonminimally cou- 
pled to gravity, containing higher order curvature invariant term. Thus, along with the lapse function N{t) , an 
additional Lagrange multiplier is required in connection with the non-dynamical Scalar field, leading to classically 
constrained gravity theory. Study of Quantum cosmological aspects of such a constrained theory may reveal more 
important cosmological implications, eg., in a recent work, such an attempt was made by Gabadadze and Shang 
[6] for a different classically constrained theory of gravity, which has been found to admit solutions, absent from 
general theory of relativity. However, the new solutions have some wonderful features like,-the spatially flat de- 
Sitter Universe can be created from nothing, has boundaries, with vanishing total energy, etc. [J. So, it appears 
worth to study the quantum cosmological aspects of the theory under consideration. Therefore, the constrained 
Hamiltonian dynamics of the theory under consideration has been studied along the line of Dirac's algorithm to 
find a constraint free primary Hamiltonian. Despite the fact that the theory under consideration contains con- 
straints in addition to the Hamiltonian constraint of the standard theory of gravity, no attempt has been made 
to find the true degrees of freedom and to find the expression for the Hamiltonian in the reduced phase space. 
Rather, it has been quantized in the same tune of Wheeler-deWitt equation, H\'^ >— 0. However, since the only 
justification for a quantized method lies in its success, so the corresponding semiclassical approximation has been 
treated here explicitly. 

The main motivation of this work is to study the connection of the quantum-cosmological equations for the theory 
under investigation with the Hamilton-Jacobi equations of the classical theory and in the process to point out 
the difficulty in finding the semiclassical time required to parameterize classical trajectories. It has been observed 
that it is nontrivial task to extract semiclassical time parameter for such a constrained system. It has been shown 
that the semiclassical time parameter extracted in the usual manner does not lead to classical constraint equa- 
tion. Further, the semiclassical time parameter may be found through speculation, but it is complex and requires 
Wick rotation. It has also been shown that in order to parameterize classical trajectories by an unique real time 
parameter, it is required to analyze Lagrangian constrained dynamics rather than Hamiltonian. 
In the following section we have written down the classical field equations corresponding to the action under con- 
sideration for isotropic and homogeneous Robertson- Walker space time. The Hamiltonian and the corresponding 
Wheeler-deWitt equation have been found in section 3, after analyzing the constraint of the theory. Semiclassical 
approximation has been made by expressing the wave function as, ^ = eTi^ , and then by expanding S in the 
powers of M instead of h. In the process, Hamilton-Jacobi equation has been found in section 4. In section 5, 
the usual procedure to extract the semiclassical time parameter has been attempted. However, it fails to fulfill 
the basic requirement of reproducing classical (Hamiltonian) constraint equation. Other naive techniques to find 
the time parameter have also been explored. Section 6 is devoted to follow an involved technique to find the 
unique semiclassical time parameter, which is the main motivation of the present paper. In section 7, a source 
term has been incorporated in the action and the same procedure has been followed. In the process, the functional 
Schrodinger equation has been found and the uniqueness of the semiclassical time parameter has been proved 
without ambiguity. The outcome of the present work has been summarized in section (8). Finally, we have added 
up an appendix in section 9, to remind the readers about the semiclassical approximation in quantum cosmology. 



2 Formalism 



The action presented by Mukhanov and Brandenberger [5], corresponding to a homogeneous and isotropic non- 
singular Universe, for which Weyl tensor vanishes, is given by, 

S - J [(1 + 0i)i? - (02 + V3<^i)y^(4i?^.i?A<''-i?2) + Vii4>i) + V2{(t>2)]V^ d^x, (1) 

where, Vi(0i) and V2{4>2) are the potentials corresponding to the pair of apparently non-dynamical scalar fields 
and (1)2 respectively. LCH may be realized from the above action (1), only for some appropriate choice of the 
potentials. Action (1) leads to Einstein's theory at small curvature, provided, Vi{(j)i) ~ (/)^ , at ^ 1, where, 
i = 1, 2. Further, to realize LCH, first requirement is to bound the curvature and the next is to obtain nonsingular 
solutions. Curvature is bounded provided, Vi((j)i) ~ 0i at |0i| 1 and de Sitter solution is obtained in the 
asymptotic region, provided, V2((/>2) constant, at \(f>2\ 3> 1. 
Now, taking Robertson- Walker line element, 

ds^ = -N{tfdt^ + a\t)[ — + r^{d9^ + sm^edcj)^)], (2) 

1 — Kr^ 

where N{t) is the lapse function, the Ricci scalar is given by, 

6 a 2 ^ N d 

and the next higher order curvature invariant term is, 

^/d d^ _n fc N d 



So the action takes the following form, 

S=^ [ [^(12</>i + 3%/l2(^2 - 6)ad2 + ^a^d<j>2 + N{12cbi + x/l202 + 6)fca - N{Vi + V2)a^]dt + Si, (3) 

oGr J iV iV 



where, the surface term is 



K = ^gij = being the trace of the extrinsic curvature. The surface term is clearly different from the one 

that appears in Einstein-Hilbert action. However, the two are the same, if (j)2 vanishes, which has got a kinetic 
term in the action, and acts as a dynamical variable. It is important to notice that the surface term is not affected 
by the presence of the scalar field (pi , since it acts as Lagrange multiplier only. In the above action N and (j)i 
act as Lagrange multipliers, resulting in a pair of constraint equations. Variation with respect to N gives the 
Hamiltonian constraint equation, but one can fix the gauge, = 1 , without loss of generality. However, we are 
not going to fix the other Lagrange multiplier (j)i , since it has been invoked to bound the curvature invariants 
for the nonsingular universe guaranteed by the field variable (f)2 ■ So, altogether we get four equations (not all 
independent), two of which are constraint equations and those obtained under the variation with respect to (fi2 
and the scale factor a are the two field equations, where we set, iV = 1 , at the end. The constraint equation 
appearing under the variation of 0i is, 

a2^a2 12' 
while the (/>2 variation equation is. 
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k V, 



k- (5) 



The Hamiltonian constraint equation, obtained by varying the action (3) with respect to N is, 



(6 - 1201 - 3\/l202)^ - VT2<}>2- + (6 + 1201 + ^1202)4 " (^i + ^2) = 0. (6) 



a ' a? 



FinaUy, varying the action with respect to the scale factor a, one gets. 



6(20i + V302-l)- + 6(20i + x/302)- + 3(20i + V302-l)^ + \/302-(60i + \/302+3)-^ + ^(Fi+F^ (7) 

In the above field equations Vl{4)i) and ¥2(4)2) denote the derivatives of the potentials with respect to 0i and 
02 respectively. For spatially flat, fc = case, the above field equations can be combined to yield a first order 
differential equation, 

"^^^ - ^"-[-7(1 - 20i)F/ + i(Fi + V2) + -4^F/02]. (8) 



The phase diagram for a particular choice of the potentials was plotted by Mukhanov and Brandenberger [2], 
showing four different classes of trajectories. In the first, the trajectory starts from the de-Sitter phase as 02 — > —00 
and evolves through to de-Sitter as 02 — > 00. In the second, for small initial values of 0i trajectory starts at 
02 = —00 reaches a turning point and returns to 02 — —00. In the third, the trajectory shows periodic solutions 
about Minkowski space-time 0i = 02 = 0. Finally, in the fourth class, trajectories starting with small 0i and 
along with 02 > 0, evolve towards de-Sitter at 02 = 00. In a nutshell, all the phase trajectories are either 
periodic about Minkowski space-time or else they asymptotically approach to de Sitter space. Hence all solutions 
are nonsingular. It is important to note that if VC vanishes, 02 turns out to be a constant and as a result V2 
vanishes and the resulting solutions again become singular. However, the condition V{' > is satisfied by all the 
potentials chosen by Mukhanov and Brandenberger . 



3 Analyzing the constraint and the Wheeler-deWitt equation 

It is important to note that in the process of developing an action that might produce nonsingular cosmological 
solutions, the determinant of the Hessian, Wij — Yl o^ dq ' corresponding to the Lagrangian of the above action 
(3) vanishes, and so the action becomes singular. Vanishing of the determinant of Hessian, signals the presence 
of constraint in the theory which should be analyzed step by step carefully. Dirac algorithm [8] is the best 
known technique to handle such constrained system and to construct the Hamiltonian. Before we proceed, let us 
recapitulate a few important artefact of Dirac's algorithm. 

1. Firstly we remember that the existence of an infinite invariance group leads to first class constraints, while 
singular Lagrangians which do not possess a local gauge invariance leads to second class constraints. 

2. The rank of Wij for a system having 2N phase space variables (that does not contain time explicitly) is N . 
For a singular Lagrangian it is R < N. under this circumstances there exists a non-degenerate R x R matrix 
Waa, such that Pi = can be solved for qa as Qa = f°'{(l,'PaA^) 1 where, a — > l,....,i?, p — > R+ \,....,N . 
In the process one finds, r = N — R primary constraints, pr = gr{q,Pa), which originate from the definition of 
momentum. 

3. Now if a function F{q,p) be defined on the primary phase space (this will be defined shortly) Fp < F, then the 
restriction on F to Fp is achieved by replacing pr by gr{q,Pa)- If F vanishes identically after this replacement, 
then it is called weakly vanishing and is denoted by F « . If in addition the gradient of F also vanishes, then it 
is called strongly vanishing and is denoted by F ~ . 

4. Primary constraints are expressed as 0r(q,p) ~ 0, and the constrained Hamiltonian, defined on the constrained 
phase space Fc , as. 

He = ^PiQi - L. 

One can write down the primary Hamiltonian defined on the primary phase space Fp , as, 
where, are the Lagrange multipliers. 

5. Next one has to check if the constraints are preserved in time, i.e., 

4>r = {0r, Hp} « {(Pr, H^} + A"{0^, 0,} EE 0. 
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If the Poisson bracket of a constraint vanishes with all other at least weakly, then it turns out to be first class 
constraint. In that case, {(f>j,Hc} ~ 0, and all these first class constraints are exhausted. The constraints 
which do not vanish modulo the constraint, leads to secondary constraints Xk ■ 

6. One needs to check the consistency condition of the secondary constraints again, i.e., if the constraints are 
preserved in time. Again for the secondary constraints, if the poisson brackets of some (say I), vanish with all 
other, then they are first class and for them {xi,Hc} » 0. Otherwise they are second class in nature and the 
Lagrange multipliers are determined in the process. There are as many undetermined Lagrange multipliers as there 
are first class primary constraints. If all the constraints are second class then the primary Hamiltonian together 
with the primary and scc;ondary constraints constitute the field equations, which may be solved in principle, leaving 
no arbitrariness in the solutions. 

7. Finally, one can try to find the true degrees of freedom and express the Hamiltonian in the reduced phase 
space Tn . However, it is a very difficult tusk, if not impossible in most of the situations. Even without finding 
the reduced Hamiltonian it is possible to quantize the theory. There is no unique prescription in this regard. Any 
quantum theory may be considered to be the correct one if it can go over to the classical counterpart through a 

suitable correspondence principk\ 

With this background knowledge, let us proceed to construct the Hamiltonian of the theory. We understand 
that in order to construct the primary Hamiltonian, one has to introduce all the constraints of the theory in the 
Hamiltonian through Lagrange multiplier. Choosing = M , M being the square of the Planck mass, one finds 
in view of the action (3) under the gauge choice iV = 1 , 

P<i>2 i _ 1 / Pa, 2^1 + ^/^<j)2 





while (pi is not invertible due to the constraint, 

Ci = p^, « 0. (9) 
Since, j^O, so the above constraint vanishes weakly. Hence the primary Hamiltonian may be expressed as, 

Hp^ = i?c + ap^i = ^PiQi - L + ap^^ 

■ ^^^^^^^pI.^ - M[(12</,i + + 6)ka - {Vi + V2)a^] + ap^, = 0, 

where, a is the Lagrange multiplier. Since, {Ci, He} docs not vanish even weakly, so Ci is a second class primary 
constraint. As there are no first class primary constraint, so there will be no undetermined Lagrange multiplier of 
the theory. Now the constraint must be preserved in time, i.e., 

Ci = {Ci,Hp,} = Di = ^^+ M{12ka - l/^V) « 0. (10) 

This is a new second class constraint as long as V/' exists, since {Ci,Ci} = {Ci,Di} = MVl'a^ . As already 
pointed out that for the existence of nonsingular solutions ¥{' must not vanish, here we again observe that V/' 
must not vanish to remove arbitrariness from the primary Hamiltonian. However no such restriction is required 
for V2 ■ Again the condition that the constraint should be preserved in time leads to, 

A = {£>!, //pj = -2^2^^, - Maa^'Vl' « 0, (11) 

modulo the constraint (10). This is not a new constraint, rather it fixes the Lagrange multiplier a = — Ij^J^tt-, 

provided VC exists. Thus both the primary and the secondary constraints are second class. Hence the primary 
Hamiltonian, being free from arbitrariness can now be expressed as, 

Tl^ 2^^ - ^P^^P^^ ] 




(1201 + V12(/.2 + 6)fca - (Vi + ^2)0-^ =0. (12) 



One can now easily check that we have obtained the correct Hamiltonian in the usual manner. The only exception 
appears to be with <p\ = = — Ma'^v" P'l>2 ' which looks like a new equation. However, this is the one, we were 



in search of, since (pi can now be inverted through this equation. In any case, it is not an independent equation, 
since taking time derivative of the first (constraint) equation (4) and using the second one (5), one can arrive at 
it. Thus the Hamiltonian, (12) being free from constraints and being able to produce all the field equations, is the 
correct one. To make a comparison, let us take Einstein-Hilbert action minimal coupled to a scalar field, 

A ^ 12M J hiad^ + k^ + j^jil<P' - Vma']dt, 
and write down the corresponding Hamiltonian, 



2 2 

^ + - 6Affca + a^V{(j)) = 0. (13) 



2AM a 2a3 

It is interesting to note that the classical field equations (4) through (7) reduce to the vacuum Einstein's equations 
at any stage of cosmic evolution as = (/>2 = 0. However, the primary Hamiltonian (12) once constructed in 
view of the constrained system under investigation, is different from (13) and does not ever reduce to the vacuum 
Einstein's equation. This fact has been manifested in the linear appearance of Pa in Hp^ . The fact that here 
Pa appears linearly in the Hamiltonian, clearly differentiates all actions (corresponding to minimal and standard 
non-minimal coupling) with the present one. 

Our next attempt will be to canonically quantize the theory under consideration. There is no standard prescription 
to quantize a classically constrained system and the only justification of a particular method of quantization lies 
in its success. The meaning of the last sentence is that one should be able to find a correspondence principle to go 
over to the classical theory under semiclassical limit. We shall write down the counterpart of the Wheeler-deWitt 
equation H\'^ >= 0, for the modified theory of gravity under consideration and show in the following sections how 
to find the notion of semiclassical time so that classical field equations may be reproduced. Now, in constructing 
the Wheeler-deWitt equation corresponding to the Hamiltonian (12), operator ordering ambiguities should be 
taken care off. Some of the operator ordering ambiguities may be removed by expressing, pa q^^{a)paq{a) , 
where, q{a) is an arbitrary function of a. Hence the first term of the Wheeler-deWitt equation corresponding to 
the Hamiltonian (12), after replacing p by —ihV, turns out to be 

1 , d^^ 1 dq a* , 



\/T2a2 dad(t>2 q da d(j)2 

Clearly, first order derivative of the wavefunction ^ with respect to a, (ie., ^) does not appear, while the same 
with respect to the scalar field (ie., appears. Further, operator ordering in p^^ and p^^ appearing due 

to the presence of the second and the third terms in (12), again introduces first derivative terms (ie., and 
) in the Wheeler-deWitt equation, corresponding to the scalars 4>i . However, from the regularity argument 

[9], it follows that can be neglected at sufficiently small values of a, ie, at sufficiently early epoch. Hence, 
the Wheeler-deWitt equation may be kept free from the first derivative terms even after removing some of the 
operator ordering ambiguities (unlike the situation encountered in standard and other nonstandard theories of 
gravity), and is expressed as, 



^2 J 1 201 H- ^/3<^!>2 - 1 2VJ, d 



M VT2a2 dad(/)2 2a^ dcpl aW{' dcpidcf) 



M 



{(1201 + ^1202 + 6)fca - {Vi + ^2)0^} 



l*>=0. 



(14) 



which is independent of the operator ordering parameter q{a) . Not all quantum states j^* > of the Wheeler 
-de Witt equation (14) are allowed, since they are constrained by the quantum analogue of the classical constraints 
(10) and (11). However, since our aim is not to find the solution of (14), so we neither, are in search of true degrees 
of freedom nor incorporate the quantum analogue of the classical constraints. Rather, for the sake of comparison, 
we write down the Wheeler-deWitt equation corresponding to Einstein's gravity with a minimally coupled scalar 
field (see Appendix), 



-(ir^ + -^ ^1^) - ^Mka^ + aV(0) 



24M^da^ a da 



I* >=0, 



where, q takes care of some of the operator ordering ambiguities. We observe that, neither , nor ^ term 
appears in the Wheeler-deWitt equation (14), which may have some deep significance, not presently known. 



4 Semiclasical approximation 



Rcparamctrization invariancc of the theory of gravity leads to the Hamihonian constraint yielding Wheeler-deWitt 
equation. Despite the fact that the theory under consideration has additional constraints, we have not made any 
attempt to find the true degrees of freedom to construct the Hamiltonian in the reduced phase space F/j . Rather, 
the corresponding quantum equation (14) has been constructed in the same tune of the standard Wheeler-deWitt 
equation J?|\E' >= 0. Now the obvious question is, "does equation (14) alone represents the quantum version of the 
classical field equations (4) through (7)"? Other way round one may ask, "is it possible to set up a correspondence 
between the Wheeler-deWitt equation (14) and the classical equations (4) through (7)"?. Since, the constraint free 
Hamiltonian represents the correct one to reproduce all the classical field equations, so the answer to the question 
raised above is positive, if under a suitable semiclassical prescription, one can find a semiclassical notion of time 
to recover the Hamilton constraint equation (6) from the Wheeler-deWitt equation (14). It is a nontrivial task 
which we shall take up in this and in the following sections. 

At energy below Planck scale, the wave function can be expressed as, ^(a, 0i,(/>2) = e^^^'^''^^''^^^ Expanding, S 
in the powers of M as, 

S = MSo + Si + M-^S2 + --- , 
and inserting it in the Wheeler-deWitt equation (14), one obtains 

ti^ 1 (i d^{MSo + Si+ M-^S2 + ■■■) 1 d{MSo + Si + M-^S2 + ■■■) d{MSo + Si+ M-^^s + •••)! 




dad(f>2 da 

td^{i\ISo + S, + M-^S2 + ---) 1 /d[MSo + S, + M-'-S2 + 



+ 



2VJ, r i d^{MSo + Si+ M-^S2 + ■■■) 1 d{MSQ + Si+ M-^S2 + ■■■) d{MSo + Si+ M-^S2 + 



-M[(12(/)i + VT2<I>2 + 6) - {Vi + V2)a^] = 0. (15) 

Let us now collect expression having same powers in M . For Einstein- Hilbert action with minimally coupled 
fields one gets to the power of an expression that states that the Hamilton-Jacobi function 5*0 depends on 
three space only (see Appendix). Here the situation is quite different in the sense that we do not get expressions 
corresponding to the order . To the order , we have. 



This is the Hamilton-Jacobi equation. It can be identified with the Hamiltonian constraint equation (6) only imder 
an appropriate choice of semiclassical time parameter. Following two sections are devoted to find the semiclassical 
time parameter. 



5 Semiclassical time parameter - standard technique 

Since in view of the Hamiltonian (12) all the velocities are now invertible and so following standard technique as 
in Einstein's gravity with minimally coupled scalar field (see Appendix, 9.2), the time parameter can be found as, 

d 1 dSp d 2MV^dSo d ( 1 dSp 20i + ^/302 - 1 . 95'o \ d 
dt ~ Vu^dhd^ ~ Vl'a^ 902 90i ^ \^a'' da ~ ^ ^ J 0^' 

The problem is that, this clioic;e of time parameter does not lead to classical constraint equation (6). This is 
because, appearing in the Hamilton-Jacobi equation (16) is not obtainable from it and as a result, remains 
arbitrary. This is the source of trouble that we encounter in the present situation, to parameterize classical 
trajectories with semiclassical time parameter following usual procedure. Thus, the standard procedure does not 
work in the classically constrained system under consideration. 
One can try to find the same by expressing the time parameter as 

^ = (b^ + + d^] ^+(-^+ + — + (u^ + + n^] — 

dt \ da d(j)\ d(j)2 J da \ da 90i 902 / 90i \ da 90i 902 / 902 ' 



where, 6, c, d, j , I, m, u, v and n are arbitrary functions of a,(j)i and (j)2. It is now possible to find d,^i and 
(j)2 , using this expression as, 



da d(f)\ d(f)2 ' 
da d(j)\ d(f)2 ' 



dSo , 9^0 

■" 

oa 



d(j)\ d(j)2 

Substituting all these expressions in equation (6) and equating the coefficients of V^o etc., between the equation 
thus formed and the Hamilton- Jacobi equation (16), one can finally arrive at the following time parameter, 




With this technique of parametrization, classical constraint equation (6) is automatically reproduced, but the 
problem associated with this time parameter is that, there still exists an arbitrariness in the form of the arbitrary 
parameter n — n(a, 01, ^2) . Thus, the time parameter is not unique, and something else should be tried. 
It is interesting to note that simply by inspection one can choose a time parameter free from such arbitrariness as, 



d_ 

dt 



1 OSq d 



1 dSo d 2Vi dSo d 



+ 



12a2 d(i)2 da da d(f)2 a^VC da d(l>2 



This time parameter reproduces equation (6), in view of Hamilton- Jacobi equation (16). However, this is purely 
intuitive on one hand and is imaginary on the other. So we must find an involved technique for this purpose. 



6 Semiclassical time parameter - an involved technique 

As already mentioned, equation (4) is an additional gravitational constraint equation of the theory under con- 
sideration, since it does not contain second order derivatives. The standard technique should be to differentiate 
equation (4) and then to compare it with equation (6). In the process, if the emerging equation is again a con- 
straint equation, then it should be entered into the Lagrangian via Lagrange multiplier. The emerging constraint 
equation in this case is, 

<^i + 2yi2^-=0. (17) 



Hence, introducing this constraint equation (17) in action (3) via a Lagrange multiplier A, we obtain 



{I2(t>i + iVvi(t>2 - 6)ad2 + y/Vla^a^2 + (12^1 + y/Vl(j)2 + 6)A:a - {Vi + ^2)0^ - A(0i + 2^12:^-) 



dt. 
(18) 



The canonical momenta are found from the action (18) as. 



V 

Pa = M[2(12(/)i + 3Vl2<^2 - 6)ad + VT2a'^(l)2 - 2VT2X^] ; p^, = -MX ; p^, = VUMa^d . (19) 

a Vi 

We don't write down the field equations since we are not going for classical solutions. However, it is important to 
note that since variation with respect to A gives back the constraint equation (17) and in view of (19) the classical 
constraint equation (6) remains unchanged, so, the primary Hamiltonian (12), the Wheeler-deWitt equation (14) 
and the Hamilton-Jacobi equation (16) are systematically reproduced. Further, identifying canonical momenta 



with corresponding derivatives of the Hamilton- Jacobi function, equation (12) can be found again from Hamilton- 
Jacobi equation (16). Now in order to obtain equation (6) from (16), let us use equations (17) and (19), which 
gives, 



■ _ _ 1 , dSo . d_ _ 1 dSo d 

. __2^p^ d_ __^fdSo._d_ 



P2 = 



_ 1 Pa _ 2(1)1 + V34>2-l p^^ _ 2y/VlV^ p^, 



s/Vl^M a3 M aV{' M 



dt^'^' ~ [Vl2a^^ da ^ ^dcj>2' aVC ^ dcj>^> j 8(^2' ^ ' 

Equations (20), (21) and (22) are now combined to yield the correct and unique semiclassical time parameter, 

d_ 1 dSo d 2V'dSo d ( 1 dSo 2cf>^ + ^/^cl>2-l dS^ 2^mVi ,^So\ _d_ 

dt ./Ua?^d(i>2'da a^i'^ dci>2' dci>^^ \^a?^ da ' a^ ^8(1)2' aV{' ^8^i'j8^2'^ ' 

One can now easily obtain and from above time parameter (23) and substituting these in equation 

(16), one can reproduce equation (6). To check whether we have found the correct and miiquc semiclassical time 
parameter let us continue by including a source term (in the form of a dynamical scalar field) in the model. 

7 Including a source term 

In this section we consider an additional source term in the action (1), in the form of a minimally coupled scalar 
field with Lagrangian density, 

Lm = -^[la„c7,0g''^ + U{a)], 

where, a is the scalar field and U{a) is an arbitrary potential. The action (1) now reads (taking M = ^ , as 
before), 

S = M J (1201 + 3\/T2<^2 - 6)^ + ^0^(1^2 + N[{12<i)^ + VU(j,2 + Q)ka - {Vi + 1/2)0^] + ^ (77 - NU{a)^ 

(24) " 

Corresponding field equations are (under variation of the above action with respect to <pi,(j)2,<J, N and a and 
setting A'' = 1), 

a2 + a2- 12' ^^^> 
a a^ k Vo 

(26) 
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a + 3-c7 + [/'(cr) =0, (27) 
a 

• 2 • '2 

M[(120i + 3\/T202 - 6)^ + 712^2- - (12(/)i + 712,^2 + 6) — + (Fi + ^2)] + ^ + U{a) = 0, (28) 
6(2.^l+^/302-l)-+6(20l+V^(/.2)-+3(2,^l+^/3</.2-l)4+V^'^2-(6</)l+\/3<^2+3)4+^^^^ 



(29) 



where, VY(0i), (02) and U'{a) are the derivatives of the potentials with respect to (l)i,(j>2 and a respectively. As 
before we differentiate equation (25) with respect to time and compare it with equation (26) to get the constraint 
equation 

.^i+2x/T2^- =0. (30) 

Next wc incorporate this equation (30) into the Lagrangian through the same Lagrange multiplier A, as before. 
In view of the Lagrangian so formed, we obtain the same set of canonical momenta (19) along with an additional 
one, corresponding to the field a, viz., 

Pa = a^&. (31) 

Variation of the Lagrangian with respect to the Lagrange multiplier A returns equation (30). So, in view of 
equations (19), (30) and (31), one can express all the velocities in terms of momenta, and thus the primary 
Hamiltonian free from arbitrariness of Lagrange multiplier is found as, 

= W^^-^^-- ''^^ " ^ ,, + ^-M[(12^.+^/l2^.+6).a-(yr+l-.)a3]+a3^(.) = 0. 

(32) 

Corresponding Wheeler- de Witt equation Hp\^ >= is. 



'M I yi2a2 dad(p2 dcg aWl' d<pid<p2 



-M- 



{(1201 + yT2<^2 + Q)ka - {Vi + ^2)0^} - + a^'Uia) ]|* >= 0. (33) 

As before, we proceed to make semiclassical approximation by expressing * [a, 01 , ^2 , cr] as = exp[^ S'(o, 0i , 02 , f)] 

expanding S in power series of M as, 5 = MSq + Si +M~^S2-\ etc. , and then substituting all these in equation 

(33). Finally, equating terms having different orders of M to zero, we find for order 

-g^=0, ie., 5o = 5o(a,0i,02). (34) 
order term leads to Hamilton-Jacobi equation for the source free nonsingular gravitational field as, 

1 _ W + Vlfe-D g, - 4|,fl§ - (12,, + + 6)*a + (V, + V.y = 0. (35) 

^/^2a^ oa 902 2a'^ 902 aW" d0i (^02 

One can easily identify this equation with the Hamilton-Jacobi equation (16). It leads to the classical source 
free gravitational constraint equation (6) under the same choice of the time parameter (23). Hence, source free 
Hamiltonian (12) can be found from it by identifying the canonical momenta with corresponding derivatives of 
So, as before. The next, ie., order term gives the following equation. 



1 r\ 



1 



d^So dSo dSi dSi dSo 



/T2a2 V ■ dad(b2 
2a3 



ih 



da d(j)2 da d 
d^So Sod Si 



2Vi 



ih- 



dSodSi _ dSi dSp 

d(f)2 dipi d4>2 dipi 84)2 



d4>l 



' d(j)2 d<f)-2 



1 f..d^Si ,dSi2 



2a3 V'^ 9a2 



+ a^U{a) ^ 0. 



(36) 



Now in order to identify this equation (36) with the functional Schrodinger equation for the matter field cr, let us 
define a function, 



/(a>i,02,cr) = D{a,(j}i,(j}2)exp{-Si). 
In view of (37), equation (36) reduces to 

1 dSo d 2Vi dSp d ( 



(37) 



1 dSo 2Vi dSo 201 + V3(f>2 - 1 dSo 



T2a2 dcj)2 da aW{ 



T2a2 da aWC 



a^U{a) 



2a3 c»ct2 

Provided, the function D(a,(j)i) satisfies the following equation, 
. 1 d'^So 2V^ 92 5*0 201 + V302 - 1 52S'o 



(38) 



1 



12a2 9a902 a^V{' 901^02 



2V^ dSp d dSn d 20i 



2a3 



I2a2^ da 902 ^84^8^' 



V302 



(39) 



Since the right hand side of equation (38) is the quantum Hamiltonian operator for the scalar field a operating on 
function / , in the background of curved space time, so it is clear that equation (38) is the functional Schrodinger 
equation for the field a , propagating in the background of curved space time, under the same choice of the time 
parameter (23). Further, under the same choice of time parameter (23), equation (39) reduces to, 



1 dD 
D~dt 



2Vi d^So ( 201 + 7302 - l\ d'^So 



12a2 dad(f>2 aW{' 901.902 



2a3 



Hence, upto this order of approximation, we have. 



vE' = lexp(^Af5o)/, 



(40) 



(41) 



where, Sq is obtainable from equation (35), / from (38) and D from (39) or (40). So at the end we observe 
that under the same choice of time parameter (23), one can parameterize the classical trajectories by reproducing 
classical source free gravitational equation (6) from the source free Hamilton- Jacobi equation (35) at one hand, 
while the M° order term (36) can be identified with the functional Schrodinger equation (38) for the field a 
propagating in the background of curved space time, on the other. This time parametrization is unique since 
equation (38) reduces to the functional Schrodinger equation only under the choice of the real time parameter 
(23) leaving no arbitrariness. This choice is correct since it fulfills both the requirements of obtaining functional 
Schrodinger equation and parametrization of the classical trajectories. 



8 Concluding remarks: 

In a series of works, Brandenberger, Mukhanov, Sornborger and others [2], [3], [4] and [5] constructed a nonsingular 
gravitational action by invoking limiting curvature hypothesis (LCH) and studied its classical aspects. The action 
contains a couple of scalar fields, one of which (02) is required to guarantee nonsingular solutions, while the 
other (01 ) binds the curvature invariants. The corresponding potentials are chosen in such a manner that for 
|0i| ^ 1, Einstein's theory is recovered at one hand and their asymptotic behaviour should be such that LCH 
may be realized on the other. In the present work, we have studied the quantum cosmological aspect of the theory 
and made a connection of the quantum cosmological equation with the Hamilton- Jacobi equation of the classical 



theory. The answer to the question that, why it is at all necessary to quantize a nonsingular Universe model, lies 
in fact that some new results have been unveiled in the process and a non-trivial technique of parameterizing the 
classical trajectories has been found. Summarily, the results are, 

1. ) Introduction of the scalar field 0i in the action (1), only introduces yet another constraint equation (other 
than the Haniiltonian constraint). Thus it is essentially a classically constrained theory of gravity and so one has 
to deal with such constraints critically, which we have done using Dirac's algorithm. 

2. ) In the absence of the scalar fields, the classical field equations reduce to those corresponding to standard theory 
of gravity. However, the Hamiltonian obtained after analyzing the constraint, does not give Einstein's theory back. 
Likewise, the Wheeler-deWitt equation is free from first derivative terms ^ , and thus is quite different from those 
obtained in view of standard and other nonstandard (nonminimally coupled) theories of gravity. 

3. ) The wave function is real for k = ±1, unlike the situation encountered in standard gravity theory (see 
appendix), where, fc = +1 leads to a complex wave function. 

4. ) The semiclassical time parameter required to parameterize classical trajectories can be found in a unique way 
only by using the technique of Lagrange undetermined multiplier in the Lagrangian constrained dynamics, rather 
than using Hamiltonian formalism. 



9 Appendix 

Paramerizing Einstein's gravity with a minimally coupled scalar field with semiclassical time: 

In the appendix wc shall recapitulate the standard technique [1^, [11], [12j, to parameterize classical trajectories 
of the Hamilton- Jacobi equation, with some appropriate choice of semiclassical time parameter, corresponding to 
Einstein's gravity with a minimally coupled scalar field, and to find the functional Schrodinger equation. 
The action for Einstein's gravity with a source term - a minimally coupled scalar field (j) : (with M = |g , which 
is different from M = ^ , used in the text) is, 

apart from a surface term. The (q) component of Einstein's equation is 
which corresponds to the Hamilton constraint equation, 



Thus, the corresponding Wheeler-deWitt equation is 

1^- 0, (44) 



2M^da^ a da dcf)'^' 2 
where, p takes care of some of the operator ordering ambiguities. 

9.1 Expansion with power series of Planck's constant h 

In the standard WKB approximation, the wave functional ^' is expressed as 4'(a,(/)) — exp [■|S'(a, ^)] , and the 
functional S{a, cf)) is expanded in the power series of the Planck's constant h as, S ^ So + hSi + h^S2 + ■ ■ ■ etc. 
Substituting all these in the Wheeler-deWitt equation (44), and collecting terms independent of h, one obtains, 

^ka^ + aV(0) = 0. (45) 

This is Einstein-Hamilton- Jacobi (EH J) equation, which is essentially equation (43), if one identifies Pa with 
and Pif, with . In order to parameterize classical trajectories, WKB time parameter is chosen in the following 
manner, 

dSo ... da I dSo . d . 1 dSo d 

— =p^ = -Maa, or, ^ = gl^^ = ' Ya^ (46) 



J- 






(dSo\ 


2M 




" 2a2 





and 



dSo 3 ; del) 1 d So . d. 1 dSo d 

— ^p,^a<t>, or, -Ql--,^^ ^e., -|, = (47) 
The classical trajectories are parameterized by a time, which is a linear combination of (46) and (47). Thus, 

dt Ma da da a^ dcj) d(t>' ^ ' 

Under this choice of time parameter (48), equation (45) produces Einstein's equation (42). This has been shown 
by Kiefer [lOj . taking, a — expa and using the unit Af = 1. It is important to note that the choice of the 
semiclassical time parameter does not involve operator ordering ambiguity. Equation (45) can in principle be 
solved for Sq and to this order of approximation one obtains ^' = exp(-|5o). Higher order terms of h produce 
fluctuations around the classical trajectory 5*0 and hence introduce corrections to the wave-functional ^' . 

9.2 Expansion with power series of Planck's mass M 

In the present work, wc have expanded S'(a, (j)) in the power series of Planck's mass, instead. The advantage is that, 
one can obtain a functional Schrodinger equation directly from it. Substituting, S — MSq + Si + M~^S2 + • • • etc., 
in the Wheeler-deWitt equation (44) and collecting terms in different orders in M , we obtain for A/^ order, 
= 0. It implies that Sq is purely a functional of gravitational field. order term gives the source free 

Einstein-Hamilton-Jacobi (EHJ) equation. 

Now, in the absence of the matter field, equations (42) and (43) reduce to 

a^ + k = 0, and, ^/a + ^ka^ = (50) 



respectively. Hence the EHJ equation (49) is essentially the vacuum Einstein equation (50), under the identification 

da 



Pa — M^M^ . To parameterize classical trajectories, we identify, 



M——=pa = -Maa, le., — = ^TT- 

oa ot a oa oa 

Further in view of equation (49), we have, 

^^=±i^/ka, ie., ^ = ±iVk-^. (52) 

Under the choice of the time parameter (52), equation (49) reduces to the vacuum Einstein's equation (50). Hence, 
expansion of S in the power series of M decouples gravity from the source and leaves it to behave classically. 
Equation (49) can be solved for , and up to this order of approximation, 

^-(0) = exp [±^Vka% (53) 

Thus the wave- functional is well behaved at the classical singularity a — > 0. The next order of approximation 
yields the following equation, 

ihfd^So,pdSo 1 d^Si\ dSodSi ^ 1 fdSi\\ ^ 

-fa V{(p) = (54) 



da"^ a da dcf)'^ J da da 2a'^ 

where, p takes care of some of the factor ordering ambiguities. Now choosing a function, /(a, (p) — D{a) exp [^] 
and substituting and from it in the above equation (54), one obtains. 



provided, D{a) satisfies the following equation, 



dSodDia) Ifd^S, ,^^^"\^(«).o. (56) 



da da 2 \ da^ a da 

Equation (55) can be identified with Tomonaga-Schwinger equation, which is essentially the functional Schrodinger 
equation for the mater field propagating in the background of curved space-time, if the semi-classical time operator 
is identified with (51). The right hand side of equation (55) represents quantum Hamilton operator Hm for the 
matter field operating on the function /(a, (f)) . Thus equation (55) reduces to 

ih^ = Ef{a, 4,) = Hmf{a, </>). (57) 
Equation (56) can be solved for D{a) as, 

D{a)^m^Qai. (58) 
Substituting S'q = ±i\/fca, in equation (56), one obtains, 

D{a) = ma^^\ (59) 
m being the constant of integration. Finally, to this order of approximation, the wave-functional takes the form, 

* = exp [^(M5o + 5i)] = exp [±^^a']f{a, 0). (60) 

Equation (55) can in principle be solved for /(a, 0) and thus explicit form of the wave- functional can be obtained. 
Now to make semiclassical approximation to the functional Schrodinger equation (55), let us express f{a,4>) = 

exp ^(a, </>)] , expand A in the power series of h as, A{a, (f)) = Ao + hA\ + fi?A2 H , substitute it in equation 

(55) and collect terms independent of h,io get. 



IdS^dAo _ J_ fdAo^ ^ 
a da da 2a^ 



a^V{4)). (61) 



Since we have already identified — ^^-^ with the semiclassical time parameter in equation (52), so the 
above equation takes the form. 



dAo 1 (dAo^"^ 



dt 2a3 V d(p 



a^V{(p) = 0. (62) 



Further, if one identifies with , it transpires that equation (62) is the Hamiltonian for the matter field in the 
background of curved space time, corresponding to the classical action, J d^x^/^Lm = 2tv^ /[^^^ ~ V{(p)]a^ dt, 
in the isotropic and homogeneous space-time under consideration. Now, since = a^4' = i so 

d_ _ 1 dAp d 
dt ~ '^'d^d^' 

and the suppressed part of the WKB time parameter appearing in equation (48) reappears. This when combined 
with equation (51), one finally obtains the semiclassical time parameter as, 

d_^_ldSod_^]_dAod_ 
dt a da da a^ d4> d(j> 
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